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ABSTRACT 



We investigate the relative importance of starbursts and AGN in nuclear activities 
of ultra-luminous infrared galaxies (ULIRGs) based on chemodynamical simulations 
combined with spectrophotometric synthesis codes. We numerically investigate both 
the gas accretion rates (m aC c) onto super massive black holes (SMBHs) and the star for- 
mation rates (m s f) in ULIRGs formed by gas-rich galaxy mergers and thereby discuss 
■ what powers ULIRGs. Our principal results, which can be tested against observations, 

| are as follows. 

c"^i . (1) ULIRGs powered by AGN can be formed by major merging between lumi- 

nous, gas-rich disk galaxies with prominent bulges containing SMBHs, owing to the 
efficient gas fuelling (m acc > 1M yr" 1 ) of the SMBHs. AGN in these ULIRGs can 
be surrounded by compact poststarburst stellar populations (e.g., A- type stars). 

(2) ULIRGs powered by starbursts with m s i ~ 100M Q yr -1 can be formed by 
merging between gas-rich disk galaxies with small bulges having the bulge-to-disk- 
ratio (/b) as small as 0.1. 

(3) The relative importance of starbursts and AGN can depend on physical prop- 
erties of merger progenitor disks, such as /b, gas mass fraction, and total masses. For 
example, more massive galaxy mergers are more likely to become AGN-dominated 
ULIRGs. 

(4) For most models, major mergers can become ULIRGs, powered either by 
starbursts or by AGN, only when the two bulges finally merge. Interacting disk galaxies 
can become ULIRGs with well separated two cores (> 20kpc) at their pericenter when 
they are very massive and have small bulges. 

(5) Irrespective of the choice of model, interacting/merging galaxies show the 
highest accretion rates onto the central SMBHs, and the resultant rapid growth of the 
SMBHs occur when their star formation rates are very high. 

Based on these results, we discuss an evolutionary link between ULIRGs, QSOs 
with poststarburst populations, and "E+A" galaxies. We also discuss spectroscopic 
properties (e.g., H/3 luminosities and line ratio of [Oiii]/H/3) in galaxy mergers with 
starbursts and AGN. 

Key words: galaxies: active- galaxies: starbursts - galaxies: nuclei - galaxies: inter- 
actions - (galaxies): quasars: general 
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1 INTRODUCTION 

The observation that ultra-luminous infrared galaxies 
(ULIRGs, defined as those with infrared luminosities greater 
than 10 12 Lq) show a mixture of two distinct types of nu- 
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clear activities, namely starbursts and active galactic nuclei 
(AGN), has led to many observational studies of their for- 
mation and evolution processes (e.g., Sanders et al. 1988; 
Solomon et al. 1992; Soifer et al. 1986; Clements et al. 
1996; Murphy et al. 1996; Sanders & Mirabel 1996; Gao 
& Solomon 1999; Trentham et al. 1999; Veilleux et al. 1999; 
Scoville et al. 2000; Surace et al. 2000; Bushouse et al. 2002; 
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Tacconi et al. 2002; Farrah et al. 2003; Armus et al. 2004; 
Imanishi & Terashima 2004; Colina et al. 2005; Iwasawa 
et al. 2005). For example, Sanders et al. (1988) proposed 
that ULIRGs formed by gas-rich galaxy mergers can finally 
evolve into QSOs after the removal of dust surrounding QSO 
black holes. Spectroscopic properties of ULIRGs have been 
extensively discussed in terms of the relative importance of 
starbursts and active galactic nuclei (AGN) in the energy 
budget of ultra-luminous infrared galaxies (e.g., Genzel et 
al. 1998; Lutz et al. 1998). 

These observations have so far raised many questions, 
the most significant being: (1) whether all ULIRGs evolve 
into QSOs, (2) what mechanisms are responsible for trig- 
gering starbursts and AGN obscured heavily by dust in 
ULIRGs, (3) what determines the relative importance of 
starbursts and AGN in spectral energy distributions (SEDs) 
of ULIRGs, (4) whether there is an evolutionary link be- 
tween starbursts and AGN in ULIRGs, and (5) whether 
there can be physical relationships between low redshift 
(\ow-z) ULIRGs and high-z dust-enshrouded starbursts and 
AGN at intermediate and high redshifts recently revealed 
by SCUBA (Submillimeter Common-User Bolometer Array) 
(e.g., Barger et al. 1998; Smail et al. 1997, 1998, 1999; Blain 
et al. 1999). Morphological studies of ULIRGs revealed that 
they show strongly disturbed morphologies indicative of vi- 
olent galaxy interaction and merging. Previous theoretical 
studies have tried to answer the above five questions in the 
context of gas fuelling to the central region of galaxy merg- 
ers (See Shlosman et al. 1990 for more general discussions 
on fuelling mechanism in galaxies). 

Physical mechanisms responsible for the formation of 
starbursts in galaxy mergers have been investigated by many 
authors (e.g., Olson & Kwan 1990; Barnes & Hernquist 1991; 
Noguchi 1991; Mihos & Hernquist 1994, 1996; Gerritsen & 
Icke 1997). For example, Olson & Kwan (1990) suggested 
that high velocity disruptive cloud-cloud collisions, which 
are more prominently enhanced in mergers, are responsible 
for the observed high star formation rates in galaxy merg- 
ers. Although these previous numerical studies provided 
some theoretical predictions on star formation rates (SFRs) 
and their dependence on the initial physical parameters of 
galaxy merging (e.g., bulge-to-disk-ratio and gas mass frac- 
tion) , they did not investigate both SFRs and accretion rates 
(ARs) onto the central super-massive black holes (SMBHs) 
simultaneously. Therefore, they did not provide useful theo- 
retical predictions on the formation and evolution of AGN, 
or on a possible evolutionary link between starbursts and 
AGN in ULIRGs. 

Physical processes of gas fuelling to the central SMBHs 
in galaxy mergers have been investigated by a number of 
authors (Bekki & Noguchi 1994; Bekki 1995; Di Matteo et 
al. 2005; Springel et al. 2005a, b). Using dynamical simu- 
lations with rather idealized modeling of gas dynamics and 
star formation, Bekki & Noguchi (1994) first investigated 
both SFRs and ARs in merging galaxies and found that 
SFRs become very high at the epoch of the coalescence of 
the cores of two merging galaxies, whereas ARs attain their 
maxima only after the coalescence. Recently, Springel et al 
(2005a) have performed more sophisticated, high-resolution 
SPH simulations including feedback effects of AGN on the 
interstellar medium (ISM), and thereby demonstrated that 
AGN feedback can be quite important for global photomet- 



ric properties of elliptical galaxies formed by major galaxy 
merging. These previous models however did not discuss the 
latest observational results of ULIRGs, partly because their 
model do not allow authors to investigate photometric and 
spectroscopic properties of dusty starbursts and AGNs in 
galaxy mergers. 

The purpose of this paper is thus to investigate simul- 
taneously both SFRs and ARs of merging galaxies in an 
self-consistent manner and thereby try to address the afore- 
mentioned questions related to the origin of ULIRGs. We 
particularly try to understand (1) physical conditions re- 
quired for galaxy mergers to evolve into ULIRG with AGN 
(or starbursts), (2) key factors which determine the relative 
importance of starbursts and AGN, and (3) epochs when 
mergers become ULIRGs with AGN. We develop a new 
model in which the physics of star formation (including gas 
consumption and supernovae feedback by star formation), 
the time evolution of accretion disks around SMBHs, and 
the growth of SMBHs via gas accretion from the accretion 
disks are included. By using this new model, we show (1) 
how SFRs and ARS in merging galaxies evolve with time, 
(2) how they depend on galactic masses, mass ratios of two 
merging spirals, and bulge-to-disk-ratios of the merger pro- 
genitor spirals, and (3) how SMBHs grow in the central re- 
gions of starbursting mergers. We also show emission line 
properties of galaxies with starbursts and AGNs by combin- 
ing the results of the simulated SFRs and ARs with spectral 
evolution codes. 

Although previous numerical simulations combined 
with spectrophotometric synthesis codes have already de- 
rived SEDs of purely starburst galaxies obscured by dust 
(Bekki et al. 1999; Bekki & Shioya 2000, 20001; Jonsson et 
al. 2005), they did not discuss at all the spectrophotomet- 
ric properties of galaxies where starbursts and AGN coexist. 
Therefore our new way of deriving spectral properties based 
on simulation results enables us to answer some key ques- 
tions raised by recent large, systematic survey of AGN (e.g., 
Kauffmann et al. 2003), such as why a significant fraction 
of high-luminosity AGN have the Balmer absorption lines. 
Previous one-zone spectroscopic models discussed what con- 
trols emission and absorption line properties of galaxies with 
starbursts and AGN (Baldwin, Phillips & Terlevich 1981; 
Veilleux & Osterbrock 1987; Kewley et al. 2001; Dopita et al. 
2006). The present simulations allow us to discuss this point 
based on the results of SFRs and ARs derived by chemody- 
namical simulations with growth of SMBHs. 

The plan of the paper is as follows: In the next sec- 
tion, we describe our numerical model for calculating SFRs 
and ARs in merging galaxies. In §3, we present the numer- 
ical results on the time evolution of SFRs and ARs and its 
dependences of model parameters. In this section, we also 
show emission line properties of galaxies mergers with star- 
bursts and AGN. We discuss the present results in terms 
of formation and evolution of ULIRGs and QSOs in §4. We 
summarise our conclusions in §5. 



2 THE MODEL 

Since the numerical methods and techniques we employ for 
modeling the chemodynamical and photometric evolution of 
galaxy mergers have already been described in detail else- 
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Table 1. Model parameters 



Model no. 


M d (x 1O 1O M ) 


k a 


h b 


ni2 C 


orbital type 


^sf ,max 


^acc.max ^ 


Comments 


Ml 


6.0 


0.2 


0.5 


1.0 


FI 


2.6 x 10° 


2.5 x 10° 


standard 


M2 


6.0 


0.2 


0.5 


1.0 


FI 


6.4 x 10 2 


x 10° 


no accretion onto SMBHs 


M3 


0.15 


0.2 


0.5 


1.0 


FI 


1.3 x lO" 1 


4.0 x 10~ 4 




M4 


3.0 


0.2 


0.5 


1.0 


FI 


5.0 x 10° 


7.0 x 10- 1 




M5 


30.0 


0.2 


0.5 


1.0 


FI 


3.0 x 10 2 


9.2 x 10 1 




M6 


0.15 


0.2 


0.1 


1.0 


FI 


2.1 x lO" 1 


6.9 x 10" 6 




M7 


3.0 


0.2 


0.1 


1.0 


FI 


5.1 x 10° 


8.2 x 10~ 3 




M8 


6.0 


0.2 


0.1 


1.0 


FI 


9.5 x 10° 


3.0 x 10~ 2 




M9 


30.0 


0.2 


0.1 


1.0 


FI 


9.8 x 10 1 


9.5 x 10" 1 




M10 


6.0 


0.2 


0.5 


1.0 


HI 


2.2 x 10 1 


2.3 x 10° 




Mil 


6.0 


0.2 


0.5 


1.0 


RR 


2.2 x 10 1 


3.6 x 10° 




M12 


6.0 


0.2 


0.5 


1.0 


BO 


3.3 x 10 1 


4.1 x 10° 




M13 


6.0 


0.2 


0.5 


0.1 


BO 


6.0 x 10° 


1.0 x 10~ 2 


LSB minor merger 


M14 


6.0 


0.2 


0.5 


0.3 


BO 


6.8 x 10° 


1.0 x lO" 1 


unequal-mass merger 


M15 


6.0 


0.2 


0.5 


0.1 


BO 


8.5 x 10° 


2.0 x lO" 1 


HSB minor merger 


M16 


6.0 


0.02 


0.5 


1.0 


FI 


6.8 x lO" 1 


4.9 x 10" 3 


gas poor 


M17 


6.0 


0.05 


0.5 


1.0 


FI 


1.8 x 10° 


1.0 x lO" 1 




M18 


6.0 


0.1 


0.5 


1.0 


FI 


8.3 x 10° 


1.6 x 10° 




M19 


6.0 


0.2 


0.02 


1.0 


FI 


1.1 x 10 2 


2.9 x 10" 3 


smaller bulge 


M20 


6.0 


0.2 


1.0 


1.0 


FI 


6.0 x 10 1 


1.6 x 10° 


bigger bulge 


M21 


6.0 


0.05 


1.0 


1.0 


FI 


3.1 x 10° 


6.1 x lO" 1 


bigger bulge, gas poor 


M22 


30.0 


0.2 


0.1 


1.0 


TI 


5.7 x 10 1 


2.4 x 10° 


tidal interaction 



a initial gas mass fraction 

^ mass ratio of bulge to disk 

c mass ratio of merging two disks 

^ maximum star formation rate (Mq yr -1 ) 

e maximum accretion rate (Mq yr~ x ) 



where (Bekki & Shioya 1998, 1999), we give only a brief 
review here. 

2.1 Progenitor disk galaxies 

The progenitor disk galaxies that take part in a merger are 
taken to have a dark halo, a bulge, and a thin exponential 
disk. Their total mass and size are Md and 7?<j, respectively. 
From now on, all masses are measured in units of Md and 
distances in units of 7?<j, unless otherwise specified. Velocity 
and time are measured in units of v = (GMd/Rd) 1 ^ 2 and tdyn 
= (Rd/GMd) 1 ^ 2 , respectively, where G is the gravitational 
constant and assumed to be 1.0 in the present study. If we 
adopt M d = 6.0 x 10 10 M and R d = 17.5 kpc as fiducial 
values, then v = 1.21 x 10 2 kms _1 and i dyn = 1.41 x 10 8 

yr- 

We adopt the density distribution of the NFW halo 
(Navarro, Frenk & White 1996) suggested from CDM simu- 
lations: 

P[r) = (r/r s )(l + r/r s y> (1) 

where r, p s , and r B are the spherical radius, the characteristic 
density of a dark halo, and the scale length of the halo, 
respectively. The dark matter distribution is truncated at 
r = 10r s corresponding to r2oo in the NFW. The value of 
r s (0.8) is chosen such that the rotation curve of the disk 
is reasonably consistent with observations. The bulge has a 
density profile with a shallow cusp (Hernquist 1990): 

p(r) oc r _1 (r + a bu i go ) -3 , (2) 

where a b uig e is the scale length of the bulge. The ratio of 



a bulge mass (Mb) to a disk mass (Md) in a disk is re- 
garded as a free parameter and represented as / b . We deter- 
mine the bulge scale length, abui gc , for a given M b based on 
the Faber-Jackson relation (Faber &: Jackson 1976) and the 
virial theorem. The bulge mass and its compactness can con- 
trol the bar formation in the disks and thus the strength of 
starbursts in mergers. The bulge contains a SMBH with the 
mass (Msmbh) following the observed relation (Magorrian 
et al. 1998); 

Msmbh = 0.006M b = 0.006/ b M d . (3) 

The model for the time evolution of Msmbh is described 
later. 

The dark matter to disk mass ratio is fixed at 9 whereas 
the bulge to disk ratio is assumed to be a free parameter 
(/ b ). The radial (7?) and vertical (Z) density profiles of the 
disk are assumed to be proportional to exp(— R/Ro) with 
scale length Ro — 0.2 and to sech 2 (Z/Zo) with scale length 
Zo = 0.04 in our units, respectively. In addition to the rota- 
tional velocity attributable to the gravitational field of the 
disk and halo components, the initial radial and azimuthal 
velocity dispersions are added to the disk component in ac- 
cordance with the epicyclic theory, and with a Toomre pa- 
rameter value of Q = 1.5 (Binney & Tremaine 1987) . The 
vertical velocity dispersion at a given radius is set to be 0.5 
times as large as the radial velocity dispersion at that point, 
as is consistent with the trend observed in the Milky Way 
(e.g., Wielen 1977). 

The ratio of Ro to Rd and that of r B to Ro are fixed at 0.2 
and 4.0, respectively, for all disk models with different Md. 
Since we adopt the scaling relation of (i s oc M d 0,5 (Kauff- 
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Figure 1. Mass distributions projected onto the rr-j/ plane for the standard model. For convenience, stellar particles (old stars) and 
gaseous ones are shown in magenta (i.e. dark matter halo particles are not shown). Big green dots represent the locations of SMBH1 and 
2. Time (T), SFRs (rh B f in units of Mqjt -1 ), ARs (m ac in units of Mqjt - 1 ) and the simulation scale are shown at upper left, upper 
right, lower right, and lower left, respectively, for each frame. Here time T represents the time that has elapsed since the simulation 
starts. Note that ARs can become very high (> IMoyr - 1 ) when the two bulges finally merge (i.e. when the two SMBHs become very 
close with each other). 



mann et al. 2003b), where fi s is the mean stellar surface 
density of a disk (described later in 2.6), Ra — C s x Md ' 25 
(or Ro oc Md ' 25 ). The normalization factor C s is deter- 
mined such that Ro = 3.5 kpc for M d = 6.0 x 10 10 M©. 
Thus the scale lengths of disks are different between models 
with different Md- 



2.2 Star formation rates 

The disk is composed both of gas and stars, with the gas 
mass fraction (/ g ) being a free parameter and the gas disk 
represented by a collection of discrete gas clouds that follow 
the observed mass-size relationship (Larson 1981). All over- 
lapping pairs of gas clouds are made to collide with the same 
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restitution coefficient of 0.5 (Hausman & Roberts 1984). The 
gas is converted solely into field stars: we do not consider the 
formation of globular clusters (GCs). Field star formation is 
modeled by converting the collisional gas particles into col- 
lisionless new stellar particles according to the algorithm of 
star formation described below. We adopt the Schmidt law 
(Schmidt 1959) with exponent 7 = 1.5 (1.0 < 7 < 2.0, Ken- 
nicutt 1998) as the controlling parameter of the rate of star 
formation. The amount of gas consumed by star formation 
for each gas particle in each time step is given by: 



Pg oc p g 



(4) 



where p g is the gas density around each gas particle. The 
coefficients in the law are taken from the work of Bekki 
(1998, 1999): The mean star formation rate in an isolated 
disk model with M d — 6.0 x 10 10 M Q and the gas mass 
fraction of 0.1 for 1 Gyr evolution is ~ 1 M© for the adopted 
coefficient (thus consistent with the observed star formation 
rate in the Galaxy; e.g., van den Bergh 2000). These field 
stars formed from gas are called "new stars" (or "young 
stars") whereas stars initially within a disk are called "old 
stars" throughout this paper. The adopted star formation 
model is similar to that with Csf = 3.5 in Bekki & Shioya 
(1998). 

Chemical enrichment through star formation during 
galaxy merging is assumed to proceed both locally and in- 
stantaneously in the present study. We assign the metallicity 
of the original gas particles to the new stellar particles and 
increase the metallicity of the neighboring gas particles. The 
total number of neighboring gas particles is taken to be N gB s, 
given by the following equation for chemical enrichment: 

AM Z = {ZiR met m s +(1.0- R mct )(1.0- ^)m s j/ mot }/7V gas .(5) 

Here, AMz represents the increase in metallicity for each 
gas particle, Zi the metallicity of the new stellar particle (or 
that of the original gas particle), ii me t the fraction of gas 
returned to the interstellar medium, m s the mass of the new 
star, and t/ me t the chemical yield. The values of -R me t and 
y me t are set to 0.3 and 0.02 respectively. 

Using numerical simulations, Thornton et al. (1998) 
demonstrated that the total amount of energy that su- 
pernovae can give to the ISM ranges from « 9 x 10 49 to 
~ 3 x 10 ergs with a typical case being m 10 ergs. This 
amount is roughly 10 % of the total amount of energy of 
Type II SN and 20 % of Type I (Thornton et al. 1998). 
They also found that most of the energy of supernovae can 
be in the form of kinetic energy within the ISM. Guided by 
these previous theoretical results, we assume that 10% of su- 
pernovae energy can be converted into kinematical energy of 
gas clouds. We adopt the Salpeter IMF with the lower mass 
cut off of 0.1 M©, the upper one of 100.0 M© and the expo- 
nent of the slope equal to —2.35 (i.e. a canonical IMF). Total 
number of supernovae at each time step can be calculated 
according to the star formation rate. The more details of the 
numerical method to give kinematical energy of supernovae 
to gas clouds are given in Bekki & Shioya (1999). 

2.3 Accretion rates onto SMBHs 

AGN activity is believed to originate from sub-parsec size 
regions at the galactic nuclei, powered by the mass accretion 
onto SMBHs through accretion disks (e.g., Rees 1984). In the 



present study, we assume that SMBHs are not surrounded by 
accretion disks in initial disks and thus the accretion disks 
are assumed to form during merging. Therefore, we need 
to model (1) formation processes of accretion disks around 
SMBHs and (2) time evolution of ARs in growing accretion 
disks in order to estimate ARs in an self-consistent manner. 
Although numerous theoretical studies have already been 
made for physical properties of static accretion disks (e.g., 
Frank, King & Raine 2002), there have been no extensive 
theoretical studies on ARs in accretion disks that are form- 
ing and growing through radial gas inflow into the central 
sub-parsec-scale region of galaxies from their outer parts. 
Furthermore only a few theoretical attempts have been made 
to elucidate the formation process of gaseous tori and accre- 
tion disks around SMBHs (e.g., Bekki 2000). 

Given this lack of theoretical detail on the evolution of 
accretion disks, we adopt the following two-fold model to 
calculate ARs. For each time step of a simulation, we first 
calculate total gas mass that can be used for the formation 
of an accretion disk around a SMBH in the central region 
of a galaxy by assuming that tidal interaction between the 
SMBH and its nearby gas clouds and the resultant destruc- 
tion of gas clouds (Bekki 2000) can be responsible for gas 
supply to the accretion disk. Then, by using a reasonable 
analytical model, we calculate the time evolution of the AR 
(m acc ) onto the SMBH for a given mass of the accretion 
disk at each time step. The details of this two-fold model 
are described as follows. 



2. 3. 1 Formation of an accretion disk around a SMBH 

The mass of an accretion disk around a SMBH is assumed 
to increase as a result of gas accretion from gas clouds being 
gravitationally trapped and destroyed by the SMBH (Bekki 
2000). We estimate the total mass of an accretion disk (M a d) 
around a SMBH based on gas densities of gas clouds within 
_R acc from the SMBH. We assume that gas clouds within 7? acc 
can be used as fuel for an accretion disk. Bekki & Noguchi 
(1994) adopted Rg as i? acc , where Rg is defined as 



R G 



GMs* 



(6) 



where Msmbh is the mass of the SMBH and a is the velocity 
dispersion (or any characteristic velocity) in the background 
components, and G is the gravitational constant. This Rg 
can be estimated to be ~ lOpc for Msmbh = 10 8 M© in a 
canonical set of galaxy parameters (Bekki & Noguchi 1994). 

Since our model for M ad evolution is based on interac- 
tion between SMBHs and gas clouds, we adopt the "Bondi" 
radius (Rb) rather than Rg as i? acc - Rb is described as 



Rb = 



2GMsmbh 

Vrcl 2 



(7) 



where v Te \ is relative velocity between the SMBH and gas. 
We assume that v IC \ is equivalent to the central velocity dis- 
persion of a bulge in each model. Therefore, i? acc is initially 
determined by the bulge mass of Mb (= /bMj) in a model 
owing to the adopted relation of Msmbh = 0. 006Mb rela- 
tion. 

Next we estimate total mass of gas clouds within i? acc 
and thereby derive a local gas density (p g ) around a SMBH. 
Guided by theoretical predictions by Hoyle & Lyttleton 
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(1941), Bondi (1952), and Ruffert & Arnett (1994), we cal- 
culate the accretion disk from gas clouds (rh c i) as follows: 



rh c \ = C B x p g , 



(8) 



where Cb is linearly proportional to Msmbh and v IC \~ A . 
Cb is therefore chosen according to Msmbh and Mb (and 
hence M d and /t,) in each model. Time evolution of M ad is 
determined by solving the following equation in terms of m c i 
and m acc : 

M ad = m c i - m acc . (9) 

7n acc represents the mass accretion rate onto a SMBH and 
we describe the way to estimate m acc later. Equation (9) 
thus means that if there is no supply of gas from outer part 
of a galaxy into the nuclear accretion disk, M a< j gradually 
decreases owing to consumption of gas within the disk. 



2.3.2 Time evolution of accretion rates 

Based on M a d, we estimate AR (rh acc ) at each time step in a 
simulation. We adopt a gas-pressure dominated standard a- 
disk with the conversion rate of accreted mass to energy(e) 
equal to 0.1 and follow the relation between M a< j and rh acc 
for a given Msmbh shown in Liu (2004): 



in- 



= 0.1(- 



M, 



ad 



5/3 

-) M yr" 



(10) 



7.7 x 1O 7 M 

We assume that rh acc should not exceed the Eddington ac- 
cretion rate, which is described as: 



"lEdd = 2.3( 



Msn 



-)M yr-\ 



(11) 



10 8 M Q 

Thus, if m aC c > w-Edd, rh acc is set to be riiEdd at every time 
step in all simulations. 

As a result of gas accretion onto a SMBH, Msmbh is 
time- dependent and its evolution is described as: 



Msmbh (t + At) = Msmbh (i) + rn acc x At, 



(12) 



where At is the time step width (corresponding to O.Olidyn) 
in a simulation. We consider that 10% of the rest mass of 
accreted gas (i.e. 0.1rh acc A£) can be converted into energy 
(-Eacc)- Although some fraction of _E acc may well be used for 
thermally and dynamically heating the ISM in galaxy merg- 
ers, it is unclear what fraction (/ agn ) of _E acc can be returned 
back to the ISM through AGN feedback effects. Accordingly, 



we compromise by assuming that / a . 



0.1, the same value 



as that derived for supernovae feedback effects (Thornton et 
al. 1998). 

The AGN feedback energy from a SMBH is assumed to 
be used for the increase of kinetic energy of gas particles 
around the SMBH. Therefore, / agn -Eacc is equivalent to the 
sum of the increase in kinematical energy of gas particles at 
each time step. The methods to give a velocity perturbation 
(directed radially away from the SMBH) to each gas parti- 
cle around the SMBH are the same as those for stellar feed- 
back effects in Bekki & Shioya (1999). The present results 
depend on / agn such that gas transfer to nuclei (thus nu- 
clear star formation and AGN fueling) can be more strongly 
suppressed in the models with larger / agn - In this paper, we 
consider that the adopted value of 0.1 is reasonable, because 
this value is similar both to that in Springel et al. (2005a) 
with / agn = 0.05 (explaining some observations) and to that 
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Figure 2. Time evolution of AR (rii acc ; top), SFR (rh s { ; middle), 
and separation of two SMBHs (Rsmbh; bottom) in the standard 
model (Ml). Note that the peak of the AR is nearly coincident 
with that of the SFR. 



by Thornton et al. (1998) for kinetic feedback effects for 
supernovae. 



2.4 Orbital configurations 

In all of the simulations of merging pairs, the orbit of the two 
disks is set to be initially in the xy plane and the distance 
between the center of mass of the two disks (r p ) is assumed 
to be ten times the disk size. The pericenter distance and the 
eccentricity are set to be the disk size and 1.0 (i.e. parabolic), 
respectively, for most of the models. The spin of each galaxy 
in a merger is specified by two angles and fa , where suffix 
i is used to identify each galaxy. 9i is the angle between the 
z axis and the vector of the angular momentum of a disk. 
fa is the azimuthal angle measured from the x axis to the 
projection of the angular momentum vector of a disk onto 
the xy plane. 

We specifically present the results of the following three 
parabolic models with different disk inclinations with re- 
spect to the orbital plane: A fiducial model represented by 
"FI" with 6>i = 0, 6 2 = 30, fa = 0, and fa = 0; a retrograde- 
retrograde model ("RR") with 61 = 180, 62 = 210, fa = 0, 
and fa = 0; and a highly inclined model ("HI") with 6>i = 60, 
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62 = 60, cj>i = 90, and fa = 0. In addition to these parabolic 
models with e p = 1, the bound orbit model ("BO") with the 
orbital eccentricity of 0.5 and with the same orbital configu- 
ration as the "FI" model is investigated. The time taken for 
the progenitor disks to completely merge and reach dynam- 
ical equilibrium is less than 16.0 in our units (~ 2.2 Gyr) for 
most of our major merger models. 

In order to compare the time evolution of SFRs and 
ARs in mergers with that of tidally interacting galaxies, we 
investigate tidal interaction models ( "TI" ) . Although we de- 
rive the results for several interaction models, we show only 
the most interesting case in the present study, since our main 
interest is on SFRs and ARs in galaxy mergers. We show a 



model with e D 



1.1 (i.e. hyperbolic), 



0, 



30, 



0, fa = and r p equal to 1.5 times the disk size. In the tidal 
interaction model, two disks do not merge at all and be- 
come separated from each other soon after their pericenter 
passage. 

All the calculations related to the above chemodynam- 
ical evolution have been carried out on the GRAPE board 
(Sugimoto et al. 1990) at the Astronomical Data Analysis 
Center (ADAC) at the National Astronomical Observatory 
of Japan. The gravitational softening parameter was fixed 
at 0.025 in our units (0.44 kpc). The time integration of the 
equation of motion was performed by using the 2nd-order 
leap-frog method. Since the masses of the bulge particle are 
set to be the same in all simulations, the initial total particle 
number in each simulation depends on the bulge mass. The 
total particle numbers for dark matter halo, bulge, stellar 
disk, and gaseous one in a model with /b = 1.0 are 60000, 
10000, 20000, 20000, respectively, in the present study. 

2.5 Emission line properties 

Our previous chemodynamical models with spectrophoto- 
metric synthesis codes for dusty starburst galaxies have al- 
ready demonstrated that major mergers between gas-rich 
spirals can become ULIRGs with L lr > 10 12 Lq, because the 
triggered nuclear starburst components can be very heavily 
obscured by dust (Bekki et al. 1999; Bekki & Shioya 2000, 
2001; Bekki et al. 2001). Since these previous studies have al- 
ready described the details of evolution from galaxy mergers 
into ULIRGs, we here do not intend to discuss the formation 
processes of ULIRGs. We instead discuss optical emission 
properties of galaxy mergers with both starbursts and AGN 
based on SFRs and ARs derived from chemodynamical sim- 
ulations. In the present paper, we discuss global, averaged 
spectral properties of galaxy mergers rather than the spatial 
difference of the properties. Two-dimensional distributions 
of emission line properties in ULIRGs will be discussed in 
our forthcoming papers (Bekki & Shioya 2005, in prepara- 
tion). 

We mainly demonstrate the time evolution of emission 
line properties of Ha, H/3, [Oiii], and [Nil] of galaxy merg- 
ers by considering the energy contribution from both ther- 
mal (i.e. starburst) and non-thermal (i.e. AGN) components. 
From the time evolution of rh s f of a merger, we first de- 
rive the SED at each time step by using stellar population 
synthesis codes. This first step is exactly the same as that 
adopted in our previous one-zone chemo-photometric galaxy 
evolution models (Shioya & Bekki 1998; 2000; Shioya et al. 
2001, 2002, 2004). Secondly, we derive the total luminosity 



600 



400 



200 







2 2 

o 

J 1 








Time (Gyr) 



Figure 3. Time evolution of the ratio (I/ a cc/i s f) of bolomctric 
luminosity coming from AGN to that from starbursts (upper) and 
bolometric luminosity from AGN (lower) in the standard model 
(Ml). Note that the time scale of the model to show the QSO- 
like luminosity (> 10 Lq) is quite short (~ O.lGyr). Note also 
that -Lacc/i'sf becomes higher (i.e. galactic nuclei dominated by 
accretion power induced activities) when L acc becomes higher. 



of the H/3 line (I/H/3) from the production rate of ionizing 
photons based on the derived SED. Thirdly, by adopting a 
typical value of the Ha-to-H/3 ratio (i.e. Hq/H/3 ~ 2.9), and 
the observed values of [Om]/H/3 and [Olll]/Ha in HII re- 
gions of nearby galaxies (Kennicutt et al. 1989), we derive 

Lua, i[oni]j an d £[nii]- 

Fourthly, we calculate the bolometric luminosity (Lboi) 
of AGN in the merger from the AR by assuming that the 
energy conversion efficiency (e) in the accretion disk around 
a SMBH is 0.1. Fifthly, by adopting a reasonable set of values 
of L[oni]/iboi = 1/300, L[ iii]/ix = 0.01, and L hol /L x = 
30 (Kraemer et al. 2004), we calculate the I/[oiii] value due 
to the AGN. Sixthly, we derive Lu_ a , L[oiii], and Lpjn]. by 
adopting typical values of Ha/H/?, [O IIl]/H/3, and [O IIl]/Ha 
in nearby galaxies with Seyfert spectra (Kennicutt et al. 
1989). Finally, we calculate the total luminosities of emission 
lines by combining the luminosities from starburst and AGN 
components. The effects of dust on spectroscopic properties 
of galaxy mergers model are not accounted for in the present 
model. 



2.6 Main points of analysis 

We mainly investigate SFRs and ARs and their dependences 
on model parameters of galaxy mergers. In order to clarify 
the relative importance of starbursts and accretion power 
induced activities in galactic nuclei, we estimate bolomet- 
ric luminosities of starbursts (L sb ) and AGN (L acc ) by us- 
ing analytic formula. L sb can be approximated as (Kennicut 
1998): 
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^ = ^u^^^&^m^ (13) 

where rj nuc is the energy conversion efficiency in nuclear fu- 
sion reactions (i.e. conversion of hydrogen to helium), and 
f ys is the fraction of massive young stars in starbursts. L acc 
can be approximated as (Shapiro & Teukolsky 1983; Frank 
et al 2002): 

where ?7 acc is conversion efficiency of rest mass energy into 
radiation in accretion disks. These equations (13) and (14) 
clearly mean that for a galaxy to show a QSO-like luminos- 
ity (~ 1O 12 L0), only small values of gas consumption rate 
(0.7M Q yr _1 ) are required if the QSO-like luminosity origi- 
nates from accretion powered activity. 

We investigate SFRs, ARs, and emission line properties 
of galaxy mergers with starbursts and AGNs, and their de- 
pendences on initial disk masses (Ma), bulge-to-disk-ratios 
(/b), gas mass fraction (/ g ), and mass ratios of two merg- 
ing disks (7112). For the models with different M d and those 
with m-2 7^ 1, we need to change masses and sizes according 
to the scaling relation of galaxies. We adopt the observed 
scaling relation by Kauffmann et al. (2003b) and derive the 
following relation: 

^«M d ' 5 , (15) 

where fi s is the mean stellar surface density of a disk. We 
determine Rd for a given Md by using the equation (15) 
and the relation of fi s oc Md/Rd 2 ■ The above scaling rela- 
tion means that less luminous galaxies show lower surface 
brightness (LSB). For convenience, the model with mi = 
0.1 (M13) including a smaller galaxy with the mass and the 
size consistent with the equation 12 is referred to as a "LSB 
minor merger". We also investigate a "HSB" minor merger 
model with U12 = 0.1 (M15) in which a smaller galaxy has a 
surface brightness 2 mag higher than the LSB minor merger 
model. 

We primarily show the results of the "standard" model 
Ml, as this model shows typical behavior for the evolution of 
SFRs and ARs. Then we show the parameter dependences 
of other models. Below, we describe the results of 22 models 
and in Table 1 summarise the model parameters for these: 
Model number (column 1), total mass of a disk (2), the gas 
mass fraction (3), the mass ratio of bulge to disk (4), the 
mass ratio wi2 of two merging disks (5), orbital types (6), the 
maximum star formation rate (7), the maximum accretion 
rate (8), and comments on the models (9). In the following 
discussion, the time T represents the time that has elapsed 
since the simulation starts. 
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Figure 4. The time evolution of SMBH1 (solid) and SMBH2 
(dotted) in the standard model. 
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Figure 5. Time evolution of projected mass density (£ ns ) of new 
stars (i.e. poststarburst populations) around SMBH1 (solid) and 
2 (dotted) of the merger (upper) and column gas density (Eg) 
around the SMBHs (lower) in the standard model. For clarity 
and comparison, only the results at nine epochs shown in Fig- 
ure 1 are described. The projected stellar density and column gas 
density are measured for particles that are located within 0.04 
in simulation units (corresponding to 0.7 kpc) around SMBHs. 
The significantly increased values of E ns after T ~ 1.8 Gyr mean 
that SMBHs are surrounded by compact poststarburst popula- 
tions formed during galaxy merging. The very large values of E g 
around T = 1.9 Gyr mean that SMBHs (thus AGN) are heavily 
obscured by metal-enriched gas (thus dust). 



3 RESULTS 

3.1 The standard model 

3.1.1 Evolution of SFRs and ARs 

Figure 1 illustrates the time evolution of morphological 
properties, SFRs, and ARs simultaneously for the galaxy 
merger in the standard model. The SFR can be moderately 
high (rh s f ~ 6.4M yr _1 ) at T — 1.86 Gyr when two disks 
can be still clearly seen as separate entities, whereas the AR 



becomes very high (m acc > lMoyr^ 1 ) at T = 1.93 - 2.00 
Gyr when two disks finally merge to form a giant elliptical 
galaxy. A disturbed outer morphology can be still seen at 
T = 2.07 and 2.56 Gyr when the AR becomes lower (i.e. 
weak AGN phases). Both the SFR and the AR become sig- 
nificantly low at T — 2.82 and 3.38 Gyr when the merger 
remnant can be morphologically identified as an elliptical 
with no peculiar fine structures (e.g., shells and plumes). 

As Figure 2 reveals, showing the time evolution of SFR 
and AR in the standard model Ml, both SFR and AR can 



be maximised when the two SMBHs become close to each 
other. This is essentially because efficient radial transfer of 
gas into the central 10 — lOOpc in the merger occurs during 
the coalescence of the two big bulges. The AR exceeds the 
0.7M©yr _1 required for QSO activity (L bo i > 10 12 L Q ) at 
T = 1.9 ~ 2.0 Gyr whereas the SFR does not exceed the 
lOOMQyr -1 required to produce a QSO luminosity. There- 
fore, this merger can be regarded as a QSO dominated by ac- 
tivities induced by accretion power of SMBHs. The epoch of 
maximum AR nearly coincides with that of maximum SFR, 
and this coincidence can be seen in most of the present mod- 
els. These results imply that (1) galaxy mergers with AGN 
activity can contain starburst components, and (2) starburst 
components in the merger would not be so easily detected 
owing to the overwhelming light from the AGN. Since mor- 
phological transformation from spirals into an elliptical is 
nearly finished at the epoch of the maximum AR in this 
model, the merger can be regarded as a forming elliptical 
with starburst and AGN components. 

As a natural result of the high ARs in the galactic nu- 
clei, the ratio of bolometric luminosity from the AGN (L acc ) 
and that from the starburst (L s b) becomes very large in 
the final phase of galaxy merging. Figure 3 shows that (1) 
L a cc/L s b becomes higher as L acc becomes higher, (2) it be- 
comes more than 100 at the epoch of maximum L acc , and 
(3) it is higher during and after the coalescence of two disk 
galaxies than before. These results suggest that the central 
starburst component in a merger can be more difficult to 
detect when the AR of the merger becomes higher. They 
also suggest that young elliptical galaxies formed by major 
galaxy merging are more likely to show spectra with AGN 
features than HII region features. We will discuss this point 
in §4.1. 

A key factor in the evolution of these systems is the 
presence of feedback from the AGN. If we construct a 
model that has gas consumption by SMBHs, but with- 
out AGN feedback, we find that the maximum rh B f and 
m acc are increased by factors of 2.8 and 113.0 respectively 
compared to model Ml (which has the feedback present). 
This model also shows a higher residual star formation rate 
(m s f ~ lOM0yr _1 ) in a sporadic way even after coalescence 
of two cores in galaxy merging. It is therefore clear from 
these results that (1) AGN feedback can suppress both (i) 
gas fueling to SMBHs and (ii) nuclear starbursts and (2) 
AGN feedback can strongly suppress residual star forma- 
tion after coalescence of two cores. Springel et al. (2005b) 
have already pointed out that AGN feedback can expel the 
remaining gas from merger remnants to shut off star for- 
mation in their sophisticated models of AGN feedback. The 
self-control of the growth of the SMBH by AGN feedback 
effects may be important for better understanding the ori- 
gin of the Magorrian relation (eg. Magorrian et al 1998). 
However, such discussion is outside the scope of this paper. 

Figure 4 shows the time evolution of Msmbh initially 
within bulges of the two merging disks. Both SMBHs grow 
quickly by a significant factor owing to efficient gas fuelling 
to the central 10 — 100 pc and the resultant formation of 
massive accretion disks around them when star formation 
rates are quite high (> lOMQyr" 1 ). The difference in the 
growth rates shown in Figure 2 is due to the fact that gas 
fuelling in the less inclined disk galaxy (i.e. galaxy 1) is more 
efficient than in the more inclined one (i.e. galaxy 2). The 
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Figure 6. B— band surface brightness ([ib) distribution at the 
epoch of the maximum AR in the standard model (Ml). 



mass of the forming elliptical is ~ 3Md corresponding to 
~ 1.8 x IO^Mq and the final combined mass of SMBH1 
and SMBH2 is 7.7 x 10 8 M Q . Therefore, the remnant el- 
liptical of this model shows MsMBH/M spn = 4.3 x 10 -3 , 
where M sp h is the total mass of the elliptical. This ratio 
of AfsMBH/Afsph is reasonably consistent with the observed 
value of 0.006 (Magorrian et al. 1998). Given the fact that 
MsMBH/M S ph = 2.0 x 10~ 3 for the model with no growth 
of SMBHs, the result shown in Figure 4 suggests that the 
growth of SMBHs during merging is quite important for el- 
liptical galaxies formed by merging to show Msmbh — M spn 
relation similar to the observed one (Magorrian et al. 1998). 

Figure 5 describes the time evolution of projected mass 
densities of new stars (Ens)) and column densities (E g ), 
which is a measure of the degree of dust extinction for a 
given metallicity (e.g., Binney & Merrifield 1998). Figure 5 
clearly shows the significantly increased values of E ns after 
T ~ 1.8 Gyr mean that SMBHs are surrounded by com- 
pact poststarburst populations formed during galaxy merg- 
ing. The very large values of E g around T = 1.9 Gyr mean 
that SMBHs (thus AGN) are heavily obscured by metal- 
enriched gas (thus dust). 

Figure 6 shows the B— band surface brightness (/xb) dis- 
tribution of the merger at the epoch of maximum AR. Two 
disks are completely destroyed to form a spheroidal compo- 
nent by violent relaxation until this epoch, and only weak 
signs of tidal disturbance can be seen in its outer stellar halo. 
Owing to the low surface brightness outer halo components 
(hb > 27 mag arcsec -1 ), this merger with a QSO-like activ- 
ity can be classified morphologically as an E if it is located 
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Figure 7. The 2D distribution of the mass fraction of new stars 
(F y ) in the standard model (Ml). The location of the SMBH1 is 
indicated by a red filled circle. Note that the SMBH is surrounded 
by young stellar populations. 

at high redshift with z > 1 (e.g., Bekki et al. 1999 for mor- 
phological properties of high z starbursting mergers). This 
result suggests that some intermediate- and high-redshifts 
QSO host galaxies with apparently spheroidal morphologies 
can be forming elliptical via dissipative major merger events. 
This result accordingly appears to be consistent with an 
observational result (Floyd et al. 2004) that QSO hosts at 
z ~ 0.4 are massive bulge-dominated galaxies. 

Figure 7 shows the two dimensional (2D) distribution 
of F y of the merger at the epoch of maximum AR, where 
F y is the mass fraction of new stars (as a proportion of 
the total number) in the central regions of the merger. This 
figure indicates that (1) the central SMBH (in the galaxy 
1) can be surrounded by young starburst or poststarburst 
components in the central 1 kpc of the merger and (2) the 
location of the SMBH is however not necessarily coincident 
exactly with the location where most of the very young stel- 
lar components are formed in the final phase of merging. 
This difference in the locations of the SMBH and the star- 
burst region does not stay significant, as the two SMBHs 
dynamically disperse the young compact starburst compo- 
nents when they become closer to each other. Most of the 
present major merger models show the coexistence of young 
starburst (or poststarburst) and AGN components in the 
central 1 kpc of mergers when morphological transforma- 
tion is nearly completed. 

We here stress that the derived coexistence of mod- 
erately strong starburst (rh s c ~ SOMeyr -1 ,) and QSO-like 
AGN frhacc ~ 3Moyr -1 ,) is due partly to gas consump- 



Figure 8. Time evolution of Lh/3 (in units of Watts) for the star- 
burst component (solid line) and the AGN component (dashed) 
in the standard model. Note that the contribution from the AGN 
strengthens significantly around T=1.9 Gyr, when the merger be- 
comes an AGN-dominated ULIRG. 

tion by the growth of accretion disks and SMBHs. Our 
model with no gas accretion onto accretion disks and SMBHs 
(model M2) shows SFR of ~ 64OM yr _1 , which is signifi- 
cantly higher than that of the standard model (See the 7th 
and 8th columns in the table 1). These comparative exper- 
iments indicate that the presence of SMBHs that can swal- 
low gas and input feedback energy can significantly influence 
SFRs in galaxy mergers. 

Owing to rapid chemical enrichment from efficient star 
formation during starburst phases of galaxy merging, the 
stellar metallicities of stellar populations that are located 
within lOOpc of SMBHs at the maximum AR (i.e. QSO 
phases) become as high as 2Zq , where Zq is the solar metal- 
licity (=0.02). The mean ages of new stars around SMBHs is 
0.85 Gyr for SMBH1 and 0.71 Gyr for SMBH2 at the maxi- 
mum AR of the merger. Given the fact that Balmer absorp- 
tion lines can become strong (thus show "E+A" spectra) ~ 1 
Gyr after dusty starbursts (e.g., Bekki et al. 2001), the above 
results strongly suggest that SMBHs in AGN-dominated 
ULIRGs can be surrounded by metal-rich and young stellar 
populations with strong Balmer absorption lines. We discuss 
an evolutionary link between ULIRGs, QSOs, and E+A's 
later in §4.3. 

3.1.2 Emission line properties 

Figure 8 shows time evolution of Lh/3 separately for the star- 
burst and AGN components. Although Lup of the starburst 
is larger than that of the AGN in the early phases of galaxy 
merging (T < 1.8 Gyr), it becomes significantly smaller than 
that of the AGN when gas fuelling to the SMBHs becomes 
efficient (T ~ 1.9 Gyr). Lh/3 of the AGN component is al- 
ways significantly larger than that of the starburst one af- 
ter the coalescence of the two bulges: owing to very minor, 




Coexistence of starbursts and AGN 11 



10' 



10" 



1 I 1 I 1 1 








V 


I.I.I 



Time (Gyr) 

Figure 9. Time evolution of the emission line ratio of [O iii]/H/3 
in the galaxy merger of the standard model. Upper and lower 
(horizontal) dotted lines represent typical values of starburst and 
AGN, respectively. Note that as galaxy merging proceeds, the 
emission line ratio evolves from starburst-like one into AGN-likc 
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Figure 10. Time evolution of the galaxy merger in the standard 
model on the [O IIl]/H/3-[N n]/Ha diagram. For clarity and com- 
parison, only the results at selected epochs are described. The re- 
sults at 5 time steps (in simulation units) are indicated by squares 
(T = 4, 8, 12, 16 correspond to 0.56, 1.13, 1.69, and 2.26 Gyr, re- 
spectively). The lower and upper crosses represent typical values 
for starbursts and AGN, respectively. Dotted and dashed lines 
represent the division between starbursts and AGN by Kewley et 
al. (2001) and Kauffmann ct al. (2003b), respectively. 
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Figure 11. Dependences of m a cc/rh s f on the initial disk masses 
(M(j) for the two sets of models with = 0.5 (solid) and /j, = 0.1 
(dotted). Both m acc and m s f are estimated at the epochs of their 
maximum values. 



S - 



r f g = 0.02 


- f g =o.05 ; 


: ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' : 




: ' ' ' ' 1 ' ' ' <■ 1 ■<-' 1 ' 1 ' : 
: f 8 =0.1 


- f,=0.2 

■ t 


|Vrr-rt~r- 






Ik-.ii .J ■ ' 



Time (Gyr) 



Time (Gyr) 



Figure 12. Time evolution of ARs (m acc ) and SFRs (m a f) for 
four models with different gas mass fraction: / g = 0.02 (upper 
left), / g = 0.05 (upper right), f e = 0.1 (lower left), and f g = 0.2 
(lower right). 



sporadic star formation after galaxy merging, L H/ 3 of the 
merger remnant is dominated by the weak AGN compo- 
nent. Figure 9 clearly demonstrates that the emission line 
ratio of [O iii]/H/3 of the merger changes from the value typi- 
cal for starbursts into that typical for AGN at T ~ 1.9 Gyr. 
These results are consistent with the result (in Figure 3) 
that L acc /Z/ S b becomes very high (> 100) at T ~ 1.9 Gyr. 

Figure 10 shows the time evolution of the galaxy merger 
on the [OlIl]/H/3-[Nll]/Ha diagram, which is often used as 
a diagnostic for determining whether spectral properties of 
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Figure 13. Time evolution of ARs (m a cc) and SFRs (m s f) for 
two models with different bulge-to-disk-ratios: /j, = 0.02 (left), 
and f h = 1.0 (right). 
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Figure 14. The same as Figure 6 but for the tidal interaction 
model (M22). 



galaxies are dominated by starbursts or AGN (e.g., Baldwin, 
Phillips & Terlevich 1981; Veilleux & Osterbrock 1987; Kew- 
ley et al. 2001). We also plot on the Figure two lines that 
demarcate the locations of AGN- and starburst-dominated 
sources. The first (dashed) line is the "extreme starburst 
classification line" from Kewley et al (2001), a theoretically- 
derived upper limit for starburst models. The second (dot- 
ted) line comes from Kauffman et al (2003b), and is based 
on a large observational set of data from the Sloan Digi- 
tal Sky Survey. Starburst spectra should reside below these 
lines, while AGN spectra should be found to the upper right. 
It is clear from Figure 10 that the merger evolves from the 
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Figure 15. The Dependence of the ratio of maximum AR 
("lace, max) to maximum SFR (m sf max ) on rh a cc,max (upper) and 
the dependence of m aC c,max on m a f max (lower) for the present 
21 models. 



middle of the starburst region toward the AGN-dominated 
region in the upper right. 

The results in Figures 8, 9 and 10 clearly show that 
there is an evolutionary link between starburst and AGN in 
terms of the spectral properties of the merger. These results 
are consistent with the derived SFRs and ARs (e.g., Figures 
2 and 3), furthermore demonstrating the capability of our 
new codes in correctly predicting spectroscopic properties of 
galaxy mergers with a coexistence of starbursts and AGN. 
Full discussions on emission and absorption line properties 
of (including spectral lines other than those discussed above, 
e.g., H7 absorption line) of galaxy mergers will be given in 
our forthcoming papers (e.g., Bekki & Shioya 2005). 



3.2 Parameter dependences 

Although the numerical results on the coexistence of star- 
bursts and AGN are similar for most of the merger models 
considered, the magnitudes of SFRs and ARs depend on 
Md, /g, fh, m 2, and orbital configuration of merging. We 
illustrate here the derived dependences and some physical 
correlations between SFRs and ARs in galaxy mergers. 



3.2.1 M d 

Figure 11 shows how the relative importance of starbursts 
and AGN in galactic nuclei of galaxy mergers depends on 
the masses of the progenitor disks. For two sets of mod- 
els with different bulge-to-disk-ratio (fh = 0.1 and 0.5), 
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rii&cc/riisf is larger for larger Mi, which means that more 
massive galaxy mergers are likely to be dominated by AGN 
rather than by starbursts. This is essentially because the in- 
terstellar gas that is radially transferred from outer parts of 
mergers can be consumed more by SMBHs than by star- 
bursts owing to the larger initial masses of SMBHs (i.e. 
Msmbh = 0.006/bAfd) in more massive mergers. We discuss 
this result later (§4) in the context of the origin of ULIRGs. 

3.2.2 / g 

Figure 12 shows that maximum SFR and ARs in galaxy 
mergers with larger gas mass fraction (/ g ) are both higher 
compared to those with smaller gas mass fraction. This is 
because the total amount of gas transfered to the vicinity of 
SMBHs is larger for the mergers with larger / g owing to a 
larger amount of gaseous dissipation in these. Figure 12 also 
shows that mergers with f g < 0.1 exhibit ARs much smaller 
than the O.7M0yr~ required for QSO activity. These re- 
sults suggest that / g is one of key parameters that deter- 
mine whether galaxy mergers show QSO activity in their 
nuclei. Irrespective of / g , the epoch of maximum SFR and 
that of maximum AR are nearly coincide with each other, 
which implies that coexistence of starbursts and AGN can 
be quite common phenomena in galaxy mergers. 

3.2.3 f b 

Figure 13 describes the time evolution of SFRs and ARs for 
the two extreme cases of mergers with fb = 0.02 (nearly 
bulge-less spiral progenitor) and fb = 1.0 (early-type M31- 
like spiral progenitor). It is clear from this figure that max- 
imum accretion rates at the epoch of maximum SFRs are 
quite different to each other in these two cases: maximum 
rh acc is 1.6 x 10M Q yr _1 for f b = 1.0 and 3.0 x 10~ 3 M Q yr _1 
for fb = 0.02. This is due to the fact that the initial Msmbh 
can determine ARs in the present model for the formation 
and the growth of accretion disks. Figure 13 also shows that 
there is little difference in the maximum SFR between the 
two models. This result suggests that fb is also a key param- 
eter which can determine whether galactic nuclei of mergers 
can be dominated by starbursts or AGN. 

The accretion radius (-Rb) within which gas clouds can 
be converted into accretion disks is initially small for small 
SMBHs in the model (M19) with / b = 0.02, owing to the 
adopted assumption of Rb oc Msmbh- Therefore, gas clouds 
have to lose a larger amount of angular momentum (with re- 
spect to SMBHs) to reach Rb during merging. As a result, 
gas is consumed by star formation rather than by the growth 
of accretion disks around SMBHs. Thus, the AR is signif- 
icantly smaller compared with other major merger models 
with bigger bulges. It should however be stressed here that 
if we relax the assumption of Rb oc Msmbh, the AR in the 
model M19 can also reach high values. 

Thus it should be stressed that it depends on the mod- 
els of accretion radius (Rb dependent on Msmbh or not) 
whether initially small MBHs (with masses of ~ 1O 6 M0 in 
the models with small fb) can grow to become SMBHs (with 
masses of ~ 1O 6 M0). Although both previous models (e.g., 
Springcl et al. 2005b) and the present one assume a rela- 
tion between Rb and Msmbh, it is not so clear (theoret- 
ically) whether there really exists such a relation around 



SMBHs in galaxies. It would be therefore safe to say that 
future very high-resolution numerical simulations, in which 
a Rb — Msmbh relation is more self-consistently determined 
from sub-pc scale gas dynamics around SMBHs, will provide 
a more robust prediction on this matter. 

3.2.4 rn 2 

Both maximum SFRs and ARs depend strongly on n%2 such 
that they are higher in mergers with larger rri2 (See the 7th 
and 8th columns of the table 1 for the models Ml, M13, 
and M14). This is because a larger amount of interstellar 
gas can be driven into the central regions of galaxy mergers 
owing to stronger tidal disturbance and the resultant larger 
amount of shock dissipation in mergers with larger 7712 . m acc 
in minor (M13) and unequal-mass (M14) mergers are well 
below IMoyr -1 , so that these merger remnants do not show 
QSO activity and thus may well be identified as low lumi- 
nosity AGN. Since these minor and unequal-mass mergers 
ultimately become SOs whereas major mergers can become 
Es (Bekki 1998), SOs are more likely to show low luminosity 
AGN activity. These results also imply that there can be a 
correlation between AGN host morphological types (e.g, Es 
or SOs) and nonthermal luminosities of AGN. 

The HSB minor merger model (M15) shows a signif- 
icantly high AR (m acc = O.2M0yr _1 ) compared with the 
LSB counterpart M13 (m aC c = O.OlMoyr -1 ), which sug- 
gests that the compactness of the smaller galaxy in a minor 
merger is a key factor for gas fuelling to the central SMBHs. 
The reason for this high AR is that the smaller galaxy of the 
HSB model can not be destroyed by the larger galaxy until 
the final coalescence of the two galaxies so that it tidally 
disturbs the ISM of the larger galaxy more strongly and for 
a longer time and thus triggers more efficient gas fuelling. 

3.2.5 Orbital configurations 

No significant dependences of time evolution of SFRs and 
ARs on orbital configurations are found for a given set of 
parameters (See the 7th and 8th columns of table 1 for the 
models, Ml, M10, Mil, and M12). Typical SFRs and ARs 
are of the order of lOMoyr" 1 and lMQyr -1 respectively, in 
these major merger models. Given the derived 7712 depen- 
dences, these results suggest major merging (mi > 0.3) is 
one of requisite conditions for QSO formation. 

3.2.6 Tidal interaction 

Both SFRs and ARs can be significantly enhanced in tidal 
interaction models, however the degree of the enhancement 
is much less remarkable compared with major merger mod- 
els. Major merger models with bigger bulges (i.e. larger fb) 
can show larger ARs, whereas tidal interaction models with 
bigger bulges do not show high ARs (~ O.7M0yr~ ). This 
is because formation of strong stellar bars, which are the 
main drivers for gas fuelling to starbursts and AGN, can 
not be formed in the bigger bulge models (e.g., fb = 0.5). 
Thus, tidal interaction models show high ARs only when 
they involve smaller bulges and larger disk masses. 

Figure 14 shows the 2D distribution of /ib in the tidal 
interaction model M22, with M d = 3.0 x 10 n M©, / b = 0.1 
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, and Msmbh = 0.012Mb at the epoch of its maximum AR 
(2.4M Q yr _1 ). Although we have investigated several rep- 
resentative tidal interaction models, only this tidal model 
with more massive SMBHs shows a sufficiently high AR 
(> O.7M0yr _1 ) to become a QSO. It is clear from Figure 14 
that interacting galaxies with QSO activity can be identified 
as two galaxies. This is quite different from major merger 
models in which a merger with a QSO activity almost al- 
ways shows a single elliptical morphology. Given the very 
limited range of parameters for tidal interaction models to 
show QSO activities, these results imply that, if they are 
formed by galaxy interaction and merging, (1) most of QSO 
hosts can be elliptical galaxies and (2) binary QSOs can be 
very rare. 

This is borne out by observational results. Mortlock et 
al. (1999) found only 16 binary QSOs, and calculated an "ac- 
tivation radius" of between 50 and 100 kpc (cf. Figure 14) . 
Assuming these were formed in a galaxy-galaxy collision, 
this implies that QSO formation occurs late in the colli- 
sion process. A more recent survey of binary QSOs from the 
Sloan Digital Sky Survey and 2dF Quasar Redshift Survey 
(Hennawi et al. 2005) found 218 quasar pairs with separa- 
tions < l/i _1 Mpc, implying a binary fraction of ~ 1 in 1000. 

3.2.7 Correlations between m s { and m acc 

Figure 15 shows that (1) there is a weak yet positive cor- 
relation between maximum SFRs (rh s f imax ) and maximum 
ARs (m acCimax ) and (2) there is a clearer correlation be- 
tween rh acCimax and m aC c,m a x/m sfimax . The above result (1) 
suggests that mergers with more pronounced starburst ac- 
tivities are likely to show more pronounced AGN ones. The 
result (2) suggests that mergers with more pronounced AGN 
activities are more likely to be dominated by AGN rather 
than starbursts. It should be however stressed here that if 
we plot m s f and rh acc from data at every time step of all 
models (including non-AGN and non-starburst phases) in 
the same way as shown in Figure 10, the derived two corre- 
lations becomes rather weak. Thus the correlations can be 
held only for mergers with strong starbursts and AGN. 



4 DISCUSSIONS 

4.1 Relative importance of starbursts and AGN 

Our numerical simulations have shown that galactic mass is 
a key factor in determining whether a forming early-type 
galaxy is dominated by starbursts or AGN. The present 
study has demonstrated that the bolometric luminosity ra- 
tio, L aC c/L a {, is larger for galaxy mergers with larger initial 
disk masses (Md). An order of magnitude estimation can 
allow us to understand this result. For global galactic star 
formation, we adopt the Schmidt law, in which rn s t oc /ig 1 ' 5 , 
where /i g is the surface gas density of a disk. We also adopted 
the Ma — /i s relation (Kauffman et al. 2003b) in the present 
study. Therefore L B f can be approximated as: 

L s { oc m s f oc /ig 1 "' oc Md°' 7j , (16) 

for /i g oc /i s in the present models with radially constant gas 
mass fraction. On the other hand, L acc can be approximated 
as: 



L acc cc m acc oc rhEdd oc Msmbh oc M b oc M d (17) 

for a given bulge-to-disk-ratio (/b). Equation 16 and 17 lead 
us to derive the following relation: 

L&cc rn&cc „, 0.25 ( 1o1 
— — oc — — oc Md . (18) 
L s( m s{ 

This relation suggests that more luminous forming early- 
type galaxies via galaxy merging are likely to be dominated 
by AGN rather than by starbursts, if merger progenitor disks 
contain a sufficient amount of gas for fuelling. 

Although the above analytically derived relation of 
"^ c ° oc Md ' 25 is qualitatively consistent with the simu- 
lations shown in in Figure 11, it is significantly shallower 
than those derived in the simulations ( m ° cc oc Md 0,95 for 
/ b = 0.1 and ^ oc M d °' 38 for f h = 0.5; See Figure 11). 
The origin of this difference might well be closely associ- 
ated with the fact that suppression of star formation from 
AGN feedback (which can enhance the relative importance 
of accretion-power-induced activity in galactic nuclei) is not 
explicitly considered in the above analytical arguments. 

It is currently less feasible to prove the above mass de- 
pendence based on the comparison between the simulation 
results and observations, because most of previous observa- 
tions focused on correlations between nuclear activities and 
the Hubble morphological types (e.g., Mouri & Taniguchi 
2004). It may well be an observationally difficult task to 
estimate separately m sf and rh acc from emission line proper- 
ties of galactic nuclei for determining L acc /L s f. We however 
suggest that future statistical studies on L acc /L s f and its de- 
pendence on galactic masses are worthwhile, because they 
can prove an example of mass-dependent evolution of galax- 
ies. 



4.2 What powers ULIRGs ? 

It has been a longstanding, remarkable problem what dom- 
inate the luminosities of ULIRGs since many observational 
studies with different wavelengths revealed possible evi- 
dences for both starbursts and AGN in ULIRGs (e.g., 
Sanders & Mirabel 1996; Lutz et al. 1998; Genzel et al. 1998; 
Tacconi et al. 2002) . Based on a mid-infrared spectroscopic 
survey of 15 ULIRGs by ISO (Infrared Space Observatory, 
Genzel et al. (1998) revealed that there is no obvious trend 
for the AGN component to dominated in the most advanced 
mergers. Lutz et al (1998) investigated the ratio of the 7.7 
fj,m PAH (polycyclic aromatic hydrocarbon) emission fea- 
ture to the local continuum for 60 ULIRGs and found that 
only about 15% of ULIRGs at luminosities below 2 x 10 12 Lq 
are powered by AGN. 

Our simulations have demonstrated that (1) initial 
galactic masses can be one of primarily important parame- 
ters that determine the relative importance of starbursts and 
AGN in galaxy mergers and (2) more massive galaxy merg- 
ers are likely to be dominated by AGN. Almost all ULIRGs 
show strongly disturbed morphological properties, which are 
the most likely to be clear signs of past major merger events 
(e.g., Sanders et al. 1988). A logical conclusion of these the- 
oretical and observational results is that if ULIRGs are more 
massive then they are likely to be dominated by AGN. Fig- 
ure 11 suggests that galaxy mergers with their progenitor 
disk masses higher than « 10 11 M Q and bigger bulges can 
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become ULIRGs with AGN. Recently Tacconi et al. (2002) 
have investigated structural and kinematical parameters of 
ULIRGs and found that ULIRGs are not so massive/bright 
as giant ellipticals. This result, combined with our simula- 
tions, suggests that ULIRGs are, on average, dominated by 
starbursts rather than AGN. This suggestion is broadly con- 
sistent with the observational results by Lutz et al. (1998) 
that about 80% of ULIRGs are found to be predominantly 
powered by starbursts. It is however not so clear why galaxy 
mergers between less luminous late-type galaxies are more 
likely to occur than those between more luminous ones at 
lower redshifts. 

4.3 Evolutionary link between ULIRG, QSOs, 
Q+A's, and E+A's ? 

Although the relative importance of starbursts and AGN in 
ULIRGs has been observationally suggested to be different 
between different ULIRGs (e.g., Genzel et al. 1998; Lutz et 
al. 1998), a significant fraction of ULIRGs have been sug- 
gested to contain starburst components (e.g., Farrah et al. 
2003). These observational results raise the following ques- 
tion: Are there any evolutionary links between ULIRGs and 
galaxies with "E+A" spectra indicative of poststarburst 
populations (e.g., A-type stars) ? This question may well 
be quite timely and important, given the fact that physical 
properties of E+A's are now being extensively investigated 
for a large number of E+A samples derived by wide field 
surveys (e.g., Blake et al. 2004; Goto et al. 2003) and by 
8m-class ground telescopes with multi-object spectrograph 
(e.g., Pracey et al. 2004). Spectral signatures of poststar- 
burst stellar populations in some QSOs (e.g., Canalizo & 
Stockton 2000; 2001) and in some ULIRGs (e.g., Poggianti 
& Wu 2000; Goto 2005) imply that there could be some close 
physical relationships between ULIRGs, QSOs, and E+A's. 
In the following discussion, QSOs with poststarburst spec- 
tra (i.e. strong Balmer absorption lines) are referred to as 
"Q+A's" just for convenience. 

The present simulations have demonstrated that 
SMBHs of galaxy mergers can be surrounded by circumnu- 
clear, compact, and young poststarburst populations when 
ARs onto SMBHs are high. This result implies that if the 
Balmer absorption lines are not significantly diluted by 
Balmer emission lines from AGNs, spectral signatures of 
poststarburst populations can be detected in galaxy merg- 
ers. Based on these numerical results, we suggest the follow- 
ing two different evolutionary paths between mergers (Mer) 
and ellipticals (Es). For SB-dominated ULIRGs that arc 
formed by merging either between less luminous disks or 
between disks with smaller bulges, strong Balmer absorp- 
tion lines can be detectable owing to less significant dilution 
of the absorption lines by emission lines from weaker AGN 
components. Therefore the evolutionary path could be; 

Mer ULIRGs E + A's^- Es. 

For AGN-dominated ULIRGs that are formed by merg- 
ing between more massive disks with prominent bulges, the 
evolutionary path could depend on whether the lifetimes of 
QSOs are shorter than the lifetimes of A-type stars. The 
present simulations have shown that the lifetime of QSOs 
(defined as the duration within which m aC c is higher than 
0.7 M Q ) is an order of ~ 0.1 Gyr (See Hopkins et al. 2005 
for possible luminosity dependent QSO lifetimes). The dilu- 



tion of the Balmer absorption lines by AGN emission can be 
significant in the very strong AGN phases so that the ab- 
sorption lines can not be detected so easily. The absorption 
lines might well be detectable either when intrinsic AGN 
luminosities become significantly smaller or when AGN are 
observed as type II (i.e. viewed from the edge of the sur- 
rounding dusty torus) . Therefore there can be the following 
evolutionary path: 

Mer ULIRGs QSOs (Q + A's) E + A' s 

Es. 

Probably, QSOs that experienced stronger starbursts in 
the gas fuelling (thus growth) processes of SMBHs are more 
likely to have detectable Balmer absorption lines. It is also 
reasonable to claim that strong Balmer absorption lines are 
more likely to be detected in type 2 Seyfert than in type 1 
owing to the less amount of dilution of stellar light by emis- 
sion from hidden broad line regions of type 2 Seyfert (with 
weaker emission lines due to dusty torus around AGN). We 
plan to investigate this point more quantitatively by nu- 
merical simulations and compare the results with already 
existing observational results (e.g., Kauffmann et al. 2003a; 
Cid Fernandes et al. 2004). 

As shown in the present chemodynamical study, the 
metallicities of stellar populations around SMBHs in galaxy 
mergers become very high (> 2Zq) owing to rapid chemical 
enrichment associated with starbursts during galaxy merg- 
ing. Furthermore, as a natural result of chemical evolution, 
the stellar metallicities around SMBHs are higher in the 
later phase of galaxy merging. These results imply that stel- 
lar populations in later AGN phases (e.g., QSOs and Liners) 
are more likely to be more metal-rich than those in ear- 
lier starburst phases for galaxy mergers with starbursts and 
AGN. Our previous chemodynamical simulations suggested 
that the abundance ratio of [Mg/Fe] after strong starburst 
of galaxy mergers can be significantly larger than the so- 
lar value due to the dominant contribution of Type II su- 
pernovae (Bekki & Shioya 1999). Therefore we suggest that 
QSOs with younger poststarburst populations can show large 
[Mg/Fe] ratios if QSOs are evolved from ULIRGs formed by 
galaxy mergers. 

4.4 Pair vs multiple mergers in ULIRGs 
formation. 

The present study has investigated mergers between two 
disk galaxies and thereby demonstrated that ULIRGs can be 
formed in the very late phase of galaxy merging when two 
galaxies nearly complete their merging. Thus the present 
model can be more relevant to ULIRGs with single cores: 
The presence of ULIRGs with multiple cores observed in 
some ULIRGs (e.g., Borne et al. 2000; Colina et al. 2001; 
Bushouse et al. 2002) can not be simply explained by the 
present pair merger models (Taniguchi & Shioya 1998). 
Previous numerical simulations showed that (1) a com- 
pact group of galaxies can be transformed into an ellipti- 
cal galaxy through multiple merging of the group member 
galaxies (Barnes 1989; Weil & Hernquest 1996), (2) repet- 
itive and multiple starbursts can be triggered by multiple 
merging of disk galaxies (Bekki 2001), and (3) the origin 
of metal-poor, hot gaseous halo of field giant ellipticals can 
be closely associated with tidal stripping of metal-poor gas 
in multiple mergers (Bekki 2001). However, previously nu- 
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merical studies did not calculate the accretion rates onto 
SMBHs and the SEDs in multiple mergers so that they could 
not provide any theoretical predictions as to (1) whether 
multiple mergers can become AGN-dominated ULIRGs or 
starburst-dominated ones, (2) in what physical conditions 
starbursts and AGN can be obscured heavily enough to be- 
come ULIRGs emitting almost all energy in infrared bands, 
and (3) what is the dynamical fate of multiple SMBHs fallen 
into the central regions of the remnants of multiple mergers. 

Thus we plan to investigate the above questions based 
on more sophisticated, higher-resolution simulations that 
allow us to study both dynamical evolution of multiple 
SMBHs and gas accretion onto the SMBHs. The results of 
these future simulations, combined with those of the present 
study, will allow us to (1) investigate what types of com- 
pact groups (e.g., spiral-rich groups) can become starburst- 
dominated ULIRGs or AGN-dominated ones (or much less 
luminous infrared galaxies) in their conversion processes into 
field elliptical galaxies via multiple galaxy merging and (2) 
discuss statistics of the observed morphological properties 
(e.g., single or multiple cores) of ULIRGs (e.g., Murphy et 
al. 1996; Zheng et al. 1999; Borne et al. 1999; 2000; Colina et 
al. 2001; Cui et al. 2001; Bushouse et al. 2002; Goto 2005). 
These simulations will also help us to understand physical 
relationships between compact group of galaxies, multiple 
mergers, ULIRGs with hot gaseous halos (Xia et al. 2002; 
Huo et al. 2004), QSOs with companion galaxies (Stockton 
& Ridgway 1991; Disney et al. 1995; Hutchings & Morris 
1995; Bahcall et al. 1997), "fossil group" with the central gi- 
ant ellipticals (e.g., Ponman et al. 1994; Jones et al. 2003). 



5 CONCLUSIONS 

We have numerically investigated both SFRs and ARs in 
forming ULIRGs via gas-rich galaxy merging in an self- 
consistent way. Dependences of the time evolution of SFRs 
and ARs on model parameters are mainly investigated. We 
summarize our principle results as follows: 

(1) ULIRGs powered by AGN can be formed by ma- 
jor merging between luminous, gas-rich disk galaxies with 
prominent bulges containing SMBHs owing to the efficient 
gas fuelling (m acc > 1M yr _1 ) to the SMBHs. AGN in 
these ULIRGs can be surrounded by compact poststarburst 
stellar populations (e.g., A- type stars). These results suggest 
that ULIRGs and QSOs can show strong Balmer absorption 
lines. 

(2) ULIRGs powered by starbursts with m s f ~ 100M Q 
yr _1 can be formed by merging between gas-rich disk galax- 
ies with small bulges having the bulge-to-disk-ratio (/b) as 
small as 0.1. As long as the accretion radii (-Rb) of SMBHs 
are proportional to the masses of the SMBHs, galaxy merg- 
ers with smaller bulges are more likely to become starburst- 
dominated ULIRGs (i.e., they can not show AGN activ- 
ity owing to a smaller amount of gas accretion onto the 
SMBHs). 

(3) The relative importance of starbursts and AGN can 
depend on physical properties of merger progenitor disks, 
such as /b, gas mass fraction, and total masses. For exam- 
ple, more massive galaxy mergers are more likely to become 
AGN-dominated ULIRGs. 

(4) For most models, major mergers can become 



ULIRGs powered either by starbursts or by AGN, when 
the two bulges finally merge. Interacting disk galaxies can 
become ULIRGs with well separated two cores (> 20kpc) 
at their pericenter only when they are very massive and 
have small bulges. These suggest that it is highly unlikely 
for interacting/merging pair of galaxies to become ULIRGs 
with double/multiple nuclei. We note, however, the results 
of Veilleux et al (2002), who found that about 7% of ULIRGs 
in their sample have nuclear separations in excess of 20kpc. 
This may suggest that ULIRGs can be formed via alternate 
routes to the major mergers examined herein. 

(5) Irrespectively of models, interacting/merging galax- 
ies show the highest accretion rates onto the central SMBHs 
and the resultant rapid growth of the SMBHs, when their 
star formation rates are very high. 

(6) ARs can become high (IMoyr™ 1 ) enough to show 
QSO-like activities (L bo i ~ 1O 12 L ) mostly in major merg- 
ers between massive disk galaxies with remarkable bulges. 
ARs however can not reach the required rates for QSOs 
(m aC c ~ O.7M0yr _1 ) in minor and unequal-mass mergers 
that form SOs. These results therefore imply that only form- 
ing elliptical via major mergers can show QSO-like activities 
whereas forming SOs (or early-type spirals with big bulges) 
via minor and unequal-mass merging show low luminosity 
AGN (e.g., type 1/2 Seyfert). 

(7) Maximum ARs (m acCjmax ) can correlate with maxi- 
mum SFRs (m s f,max) in the sense that galaxy mergers with 
higher rh sfjmax are likely to show higher rh acc , max . This sug- 
gests that mergers and ULIRGs with more pronounced AGN 
activities are likely to show stronger starburst components 
in their nuclei. The correlations can be discussed in the con- 
text of recent observational results (e.g., Goto 2005) on cor- 
relations between infrared luminosities of ULIRGs, star for- 
mation rates, and AGN luminosities (measured from [OIII] 
emission lines). 

(8) The ratio of rhacc.mEtx to m acc ,sf can correlate with 
?Ti a cc, ma x in the sense that galaxy mergers with higher 
m acc , max are likely to show higher rh acc , max /m sfimax . This 
implies that merger and ULIRGs with higher AGN (thus to- 
tal) luminosities are likely to be dominated by AGN rather 
than by starbursts. This result can be also consistent with 
recent results on AGN fraction as a function of infrared lu- 
minosities of galaxies (e.g., Goto 2005). 

(9) There could be evolutionary links between ULIRGs, 
Q+A's, QSOs, and E+A's. Galaxy mergers between less 
massive disks are more likely to evolve from starburst- 
dominated ULIRGs into E+As without experiencing QSO 
phases, whereas those between more massive disks with 
prominent bulges can evolve from AGN-dominated ULIRGs, 
to QSOs (and/or Q+A's), and finally to E+A's, if the life- 
times of QSOs are as short as ~ 0.1 Gyr. Removal of gas 
reservoir for star formation via supernovae and AGN feed- 
back could be essentially important for the above evolution- 
ary links. 

(10) Time evolution of emission line properties of galax- 
ies with starbursts and AGNs is investigated based on SFRs 
and ARs derived from chemodynamical simulations. For ex- 
ample, simulated mergers are demonstrated to evolve from 
those with smaller [O iii]/H/3 (starburst-dominated) to those 
with larger [O m]/H/3 (AGB-dominated) . It is suggested that 
strong Balmer absorption lines are more likely to be detected 
in type 2 Seyfert than in type 1 owing to the less amount of 
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dilution of stellar light by emission from hidden broad line 
regions of type 2 Seyfert. Direct comparison between the 
predicted spectrophotometric properties of galaxy mergers 
with dusty starbursts and AGNs and the corresponding ob- 
servations will be done in our forthcoming papers. 
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